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We consider a quantum system linearly coupled to a reservoir of harmonic oscillators. For finite 
coupling strengths, the stationary distribution of the damped system is not of the Gibbs form, in 
contrast to standard thermodynamics. With the help of the quantum Hamiltonian of mean force, we 
quantify this deviation exactly for a harmonic oscillator and provide approximations in the limit of 
high and low temperatures, and weak and strong couplings. Moreover, in the semiclassical regime, 
we use the quantum Smoluchowski equation to obtain results valid for any potential. We, finally, 
give a physical interpretation of the deviation in terms of the initial system-reservoir coupling. 

PACS numbers: 03.67.-a, 05.30.-d 
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I. INTRODUCTION 

Thermodynamics successfully describes equilibrium 
properties of macroscopic objects coupled to a heat bath. 
In the standard theory, the system-bath interaction is 
assumed to be small and is therefore neglected. This is 
a well-justified assumption when the coupling energy is 
small compared to the thermal energy of the system pQ. 
However, at low temperatures, quantum effects come into 
play and the interaction energy cannot be discarded any- 
more [2] • This is particularly relevant for solid-state nan- 
odevices, such as nanomcchanical oscillators, which may 
be strongly coupled to their environment [3H6]. In the 
limit of vanishingly small coupling, the damped quan- 
tum system asymptotically relaxes to the correct ther- 
mal Gibbs state [9]. However, for any finite interaction 
strength, the stationary state of the quantum system de- 
viates from a Gibbs state, in contrast to the prediction 
of thermodynamics. It has, for example, been shown 
that, at zero temperature, the coupled oscillator is in an 
excited, mixed state and not in its pure ground state jTDJ- 
rT2j . In addition, it is worth noticing that the deviation 
from a Gibbs state leads to the interesting possibility of 
entangling two noninteracting coherent states via the in- 
direct coupling to a common reservoir |13j . 

A useful quantity to describe classical systems that are 
strongly coupled to their environment is the so-called po- 
tential of mean force jUj . It was first introduced by Kirk- 
wood in his study of the average structure of liquids |15j . 
It has since proved to be an invaluable tool in chemistry 
[16| . where it is commonly employed in the investigation 
of implicit solvent models |17) and of protein-ligand bind- 
ing |18j . to name a few examples. A quantum extension 
in the form of the Hamiltonian of mean force (HMF), see 
Eq. ([5]) below, has recently been proposed in Ref. [15] . 
Potential and Hamiltonian of mean force are instrumen- 
tal in characterizing the deviations from standard ther- 
modynamics that occur in the strong coupling regime. 
However, they are difficult to evaluate explicitly in the 
general case. Here, we use the quantum HMF to quantify 
the deviation from the thermal Gibbs state for a damped 
quantum harmonic oscillator, for which we obtain exact 
analytical results. We further derive approximate expres- 



sions for arbitrary systems in the strong damping regime. 
We additionally relate the HMF to the thermodynamic 
change of energy and entropy associated with the initial 
coupling between system and bath. 

The present article is organized as follows. In Section 
[TTJ we introduce the concept of the quantum HMF and 
show how it can be used to characterize the deviation 
from a thermal Gibbs state. In Section \LU\ we determine 
the exact expression for this deviation for the case of the 
analytically solvable damped quantum harmonic oscilla- 
tor. In Section IV we provide useful approximations in 
the limit of high and low temperatures, weak and strong 
coupling. Furthermore in the semiclassical limit, we em- 
ploy the quantum Smoluchowski equation to evaluate the 
deviation from the Gibbs state for any damped quantum 
system. In Section |VI[ we finally give a physical inter- 
pretation of the HMF by establishing a relationship with 
the initial system-bath coupling. 



II. HAMILTONIAN OF MEAN FORCE 



We consider the total Hamiltonian, 



H = H S + H E 



H, 



SB 



(1) 



which describes a quantum system S coupled to a heat 
bath B via the interaction SB. The equilibrium state of 
the total system is given by the Gibbs state, 



P : 



exp(-/3g) 
Z 



(2) 



where Z — tr{exp(— /3H)} is the total partition function 
and /3 = l/(kT) the inverse temperature of the bath. 
When the interaction energy is negligible, H ~ Hs + Hb, 
the reduced density operator of the system, ps — trs{/o}, 
is also of the Gibbs form, 



PS 



exp(-/3g s ) 
Z s 



(3) 



with Zs = tis{exp(—(3Hs)}- However, for finite cou- 
pling, the interaction term cannot be neglected, and the 



2 



reduced density operator is no longer of the Gibbs form. 
It can generically be written as, 



PS 



exp(-L3H* s ) 



(4) 



where the quantum HMF is defined as the effective 
Hamiltonian \W\, 



= _ 1 tr B {exp(-/3g)} 
5 n tr B {cM-PH B )} 



(5) 



The corresponding partition function is given by Z* = 
tr s {exp(— )3Hg)}. The HMF reduces to the Hamiltonian 
of the system, H$ ~ Hs, in the limit of vanishing cou- 
pling. It is thus convenient to introduce the difference, 



AH S = m-H s • 



(6) 



The latter quantifies the deviation from a thermal Gibbs 
state. In the following, we evaluate the deviation AH$ 
explicitly for the damped quantum harmonic oscillator. 



Due to the linearity of the model, the dynamics of the 
reduced quantum system is exactly solvable. Its station- 
ary state is conveniently described in phase space using 
the Wigner function [2Tj . 



W(q,p,t) = 



2-rrh 



exp 



(12) 



where ps(x,x',t) = (x\ps{t)\x') is the position represen- 
tation of the reduced density operator. The Wigner func- 
tion of a damped harmonic oscillator is found to satisfy 
the exact quantum master equation |22j . 

g-W(q,p,t) = -^—W+—[( 1P + Muj 2 q)W} 



dp 



if- 



+D pp (t) — W + D qp lr,^\l. (i.T) 



dqdp 

The detailed expressions for the diffusion coefficients D 



!>!> 



and D qp can be found in Ref. [22] . The reduced station- 
ary phase space distribution for the Ohmic model is given 
by the Gaussian Wigner function, 



III. MICROSCOPIC SYSTEM-RESERVOIR 
MODEL 



The standard model for a damped quantum harmonic 
oscillator is given by Hamiltonian ([I]) [20], where 



He 



P 

2M 



M 



2 2 

U! q 



(7) 



is the Hamiltonian of the harmonic oscillator with mass 
M and frequency u>, 



N 

£ 

3=1 



2m 



mjiv 2 



(8) 



describes the bath of N harmonic oscillators, and 



H 



SB 



N 
.7=1 



-CjqXj 



(9) 



is the coupling term linear in the position of the system. 
The frequencies of the reservoir modes are taken to be 
equidistant, that is LUj = jA. The coupling constants Cj 
are chosen to obey a Drude-Ullersma spectrum (20) . 



c 3 = 



l2ymjMujjA 



J D 



(10) 



J D 



with damping coefficient 7 and cutoff frequency u>d. In 
the continuous limit N —> 00 (A — > 0), the bath is char- 
acterized by the Ohmic spectral density function |20j . 



7T \ ■> 



C- 71/Mwf, 
77/ — — o(u — uij) = ^ if 



(11) 



W(q,p) 



1 



2n^W)W) 



exp 



P 



2(q 2 



2 (pi 



(14) 



with exact position and momentum dispersions 



3 
3=1 



(A i -w D )V(l + ^ 



2tt 



(A 



i+i 



A,-)(A,-_i - Aj) 



1 



(P 2 



M(3u 2 ' 



(15) 



£ 

i=i 



A,V( 



1 



2tt 



(A. 



i+i 



Aj-)(Aj--i - Aj) 



Af 2 W V)(r7,M) 



(16) 



The parameters Ai are the characteristic frequencies of 
the damped oscillator and ip is the digamma function. 
In the limit of large cutoff frequencies, ljd 3> 7,w, the 
characteristic frequencies can be approximated by 23 , 



Al,2 

A3 



1± T 
2 V 4 

uj d - 7 . 



(17) 



A direct consequence of Eq. ( 16 ) is that equipartition is 
in general violated, Muj 2 (q 2 ) <(p 2 )/M, and ps is not of 
the Gibbs form. This can be understood by noting that 
due to the form of the coupling Eq. (|9| , the bath can be 
seen as continuously measuring the position of the system 
|20) . As a result, the variance of the position is reduced. 
The momentum dispersion then increases following the 
Heisenberg uncertainty relation. For a vanishingly cou- 
pled oscillator, the variances reduce to 



Moj 2 (q 2 ) 1=0 



(p 2 ) 7 =o _ fa 



M 



= — coth 
2 



(18) 



3 



Equipartition therefore holds and the reduced station- 
ary state is Gibbsian. The position representation of the 
reduced density operator is further given by 20j, 



ps(q,q') 



1 



exp 



{q + q'Y {q-q') 



l\2 



8(q 2 ) 2h 2 /{p 2 ) 



(19) 



IV. 



HMF FOR DAMPED QUANTUM 
OSCILLATOR 



The exact expression for the quantum HMF ^ for the 
damped oscillator can be obtained in the following way 
[23"] . Since the stationary state Eq. (14) is a Gaussian, 



it can be regarded as the equilibrium Wigner function of 
an effective harmonic oscillator with Hamiltonian, 



m = 



2M* 



1 ,,r* *2 2 

-M uj q 



(20) 



with mass M* and frequency uj* . The latter reads 



W(q,p) 



tanh(/3Sw*/2) 



exp 



2H* Q 



huj* coth (/3Huj*/2) 



(21) 



By comparing the two Wigner functions Eq. (|14j) and 
Eq. (21), one can identify the effective mass, 




(p 2 ) Ph 



M* 



and the effective frequency 
- 2 



1 



2 arcoth (2d) ' 



W = — arcoth (2v) 



(22) 



(23) 



Here v = y/ (p 2 ) (q 2 ) /h is the phase space volume. In 
order to estimate the effect of the finite coupling on the 
reduced state of the system, we now determine the devi- 
ation AH S given by Eq. Combining Eqs. pOj ), ( [22] ) 
and (23) for the quantum HMF, we arrive at 



AH; 



s 



2 2 

At+B q 



2 2 ' 

where we have introduced the two coefficients 



A 



B = 



'(a 2 ) 2 1 

-^rf — arcoth(2i;) 

(p 2 ) f3h y ' M 



(q 2 ) (3h 



arcoth(2w) - Muj 2 



(24) 

(25) 
(26) 



Equation (24) gives the exact expression for deviation 



AHs for the damped quantum oscillator. The coeffi- 
cients A and B are divergent for 2v = 1, since arcoth a; 
is not defined at x = 1. This happens for an undamped 



MA 



MA 




FIG. 1: (Color online) (a) Coefficient MA, Eq. (25 1 



function of the dimensionless parameters kT/(hw) and 7/w. 
A increases monotonically to zero with temperature and in- 
verse coupling strength. The deviation from a Gibbs state 
in momentum is thus maximal in the strong-coupling, low- 
temperature limit, (b) Cross sections of the (T, 7)-curve for 
kT/(hoj) = 0.5 (red lower line), 3 (blue line), 15 (green up- 
per line). In the limit of high damping, A coincides with the 
semiclassical expression given by the quantum Smol uch owski 
equation (black dashed line for kT/(hui) — 3), Eq. (371. Pa- 
rameters are ujd = 1000, uj = 1, M = 1. 



harmonic oscillator (7 = 0) in its ground state (T = 0). 
However, in this case A = B = 0, as the system relaxes 
to a Gibbs state. 

Figure [l] shows the coefficient A, Eq. ( 25 ) , as a func- 
tion of temperature and coupling strength. A is always 
negative, since — 1 < MA < 0. We observe that A de- 
creases monotonically with the coupling and the inverse 
temperature. The maximum deviation from the Gibbs 
state is thus achieved in the low-temperature, strong 
coupling limit. Moreover, for weak coupling, A van- 
ishes rapidly with growing temperature, whereas it de- 
cays much more slowly in the limit of strong coupling, 
indicating that the deviation from a Gibbs state persists 
at much higher temperatures in the latter. The coeffi- 



B 

Moj 2 




FIG. 2: (Color online) (a) Coefficient B/(Muj 2 ), Eq. (26 1 



as a function of the dimensionless parameters kT/(hui) and 
y/u). B increases monotonically to zero with temperature. 
However, it exhibits a minimum as a function of the cou- 
pling strength, (b) Cross sections of the (T, 7)-curve for 
kT/{huj) — 0.1 (red lower line), 0.5 (blue line), 1 (green 
upper line). For increasing temperatures, the dip gets less 
pronounced and eventually disappears. In the limit of high- 
damping, B coincides with the semiclassical expression given 
by the quantum Smoluchowski equation (black dashed line 
for kT/(htu) — 0.1), Eq. (451. Same parameters as in Fig. [I] 



FIG. 3: (Color online) Function arcoth(2w) (red continuous) 
and J (p 2 ) I (q 2 ) (blue dashed) appearing in the definition of 
the coefficient B, Eq. \2Q\ , as a function of the dimensionless 
parameter 7/aj. Their different behaviors lead to a minimum 
for their product (green dotted line). Same parameters as in 
Fig. [l] and kT/fiw = 0.1. 



ing coupling constants. Its product with the increasing 
square root in front thus displays a minimum (see Fig.|3|. 
Figure 2(b) additionally indicates that the minimum gets 



less pronounced for higher temperatures and eventually 
disappears in the classical limit. Finally, in Fig. |4j we 
show the average relative deviation (AHs) / (Hs), which 
vanishes in the high-temperature limit. It is worth notic- 
ing that this behavior is peculiar to the linear coupling 
model that we consider here. For the case of nonlinear 
system-bath coupling, the deviation exists even in the 
classical limit |25j . Since A(p 2 ) is in general much larger 
than B(q 2 ) (see below), the effect we have observed for B 
is suppressed, and the deviation decreases monotonically 
with increasing coupling strength. In the following, we 
derive analytic expressions for the coefficients A and B 
in various limits of interest, including the limits of low 
and high temperatures, weak and strong couplings. In 
all cases, we assume the cutoff frequency wrj to be large. 



A. High-temperature limit 



In the high-temperature limit, hu> <C kT, the posi- 



cient B, Eq. (26), is plotted in Fig. 2(a) as a function 



of temperature and coupling strength. It is also negative 
with —1 < B/(Mu> 2 ) < 0. While B shows a similar tem- 
perature dependence as the coefficient A, its dependence 
on the coupling strength differs significantly. The three 
lines highlighted in Fig. 2(a) and reproduced for clarity 
in Fig. |2(b)| show that B possesses a minimum. Hence, 
the deviation of the Gibbs state in position increases up 
to a certain value of 7 before decreasing even though the 
coupling becomes larger. The presence of a minimum 
can be understood by noting that (q 2 ) decreases with in- 
creasing coupling constant, while (p 2 ) increases [20]. As a 
result, the ratio (p 2 ) / (q 2 ) increases faster than the prod- 
uct (p 2 )(q 2 ). The hyperbolic arccotangent of the phase 
space volume in Eq. (26), therefore, decreases for grow- 



tion and momentum quadratures (15) and (16) can be 
expanded to lowest order to yield |23j . 



<<7 2 > 
(P 2 ) 



kT 
Muj 2 



MkT 



1 



1 (huj 
12 \kf 

1 h 2 {u 2 + juj D ) 

12 (kT) 2 



(27) 
(28) 



where we have made use of the series expansion of the 
digamma function Eq. ( A2 1 . We note that the coupling 



constant 7 only appears in second order in (tkj/kT). The 
two variances Eqs. (27) and (28), hence, become damp- 
ing independent in the classical limit, and equipartition 
holds, (p 2 )/M — Muj 2 (q 2 } = kT. For high temperatures, 
the phase space volume v is large and the arcoth in the 



5 




(AH S ) 
(Hs) 



FIG. 4: (Color online) Relative mean deviation (AHs) / (Hs) , 
Eq. (24 1, showing the departure from a Gibbs state as a 



function of dimensionless temperature and coupling strength. 
The maximum deviation is observed in the low-temperature, 
strong-coupling regime. Same parameters as in Fig. [I] 



coefficients A and B, Eqs. (25) and (261, can be approx- 



imated with the help of Eq. (A3). We obtain, 



A 



D 



1 

12M 

Muj 2 
12 



JUJ D ) 



(kT) 2 

2 



kf ) 



A jT 2 



(29) 

(30) 
, see 



Both coefficients approach zero quickly as 
Figs. [T] and 2(a) However, unlike B, the coefficient A 
depends explicitly on 7 and ujr>. This feature follows 
from the Ohmic nature of the damping, as seen from 
Eq. (28) for (p 2 ) . In the limit of high cutoff frequency, A 
is much larger than B. Combining Eqs. |27"[)-(30 ) , we find 
that the total deviation AHs ~ —1/T in the classical 
limit, see Fig. |4j In this regime, the stationary state 
of the system reduces to a classical Gibbs distribution, 
W(q,p) = j\fexp[-(Muj 2 q 2 + p 2 /M)/(2kT)}. 



B. Low-temperature limit 

For low temperatures, kT <C Hlu, we can approximate 
the digamma function in Eqs. 
help of Eq. 



(All 



(151 and (16) with the 



To proceed further, it is important 
to distinguish between weak and strong coupling limit. 
For weak coupling, kT <C fvy <C huj, an expansion up to 
second order in kT/huj and first order in 7/w yields |23) . 



<9 2 



(P 2 



2Mw 2 



1 - 



7 



1 



7 




2 


CO 


K7T 


7T 




(31) 



(32) 



For 7 = 0, Eqs. (31 ) and (|32| lead to the correct ground 



(H s ) = (p 2 )/2M + Muj 2 (q 2 )/2 = hu/2. At the same 
time, the phase space volume can be approximated by 
v = 1/2(1+C), where C = j(ln(ui D /uj)-l)/(uJ7r). There- 
fore, making use of Eq. ( |A4[ ), we find the coefficients, 

' kT 

' -(l-y)ln(C/2) + l 



.4 



1 

M 



B 



-Mlj 2 



'kT 

huj 



(1 + Y) ln(C/2) + 1 



(33) 
(34) 



with Y = 7ln(aj£)/w)/(w7r). Here, in contrast to the 
high-temperature limit, both A and B increase with T. 
As a result, the deviation from the Gibbs state decre ases 
with increasing temperature, see Figs. [I] and 2(a) At 
zero temperature, the coefficients A and B attain their 
minimum values, A = —1/M and B = —Muj 2 , that cor- 
respond to the maximum deviation from the Gibbs state. 

In the opposite limit of strong damping, kT <C huj <C 
hj, a lowest order expansion of the variances, Eqs. ( 15 ) 



and ( 16 ) leads to 



(<1 ) = — 77- In - 4 
irM 7 uj 

M^7 ln c^ 

7T 7 



nhrf ( kT \ ' 
3Moj 2 \hu~j~) 



(35) 



(36) 



In this regime, the phase space volume v is large and the 
arcoth in A and B, Eqs. (25) and ([261, can therefore be 



approximated by Eq. (A3). We then obtain in first order, 

1 f kT UJ 7T 



A — z — \ (37) 

M hu 7 \xv(uj£)/^i) 

B = M ^\l¥L ZL il . (38) 
[oj Huj 21n(7/w) J v ' 

Again the two coefficients increase linearly with T, see 
Figs. [T] and [2(a) 



SEMICLASSICAL LIMIT 



In the overdamped regime, the dynamics of the quan- 
tum particle simplifies considerably. In the limit 7/cj 2 3> 
(hj3, 1/7), the off-diagonal matrix elements of the reduced 
density operator of the system in the coordinate repre- 
sentation are strongly suppressed, and a semiclassical de- 
scription becomes possible. On the coarsed-grained time 
scale, t 3> I/7, the diffusion coefficients in the master 



equation (13) simplify to [2"6l |2"T] . 

= V"(q)(q 2 )- 
l(p 2 ) ■ 



D, 



qp 



(e!> 

M 



(39) 
(40) 

The position and momentum variances for the harmonic 
oscillator are given in this limit by [28], 

(9 2 > = 75177^ +A, (41) 



PMuj 2 



state mean energy of an undamped harmonic oscillator, 



&M7 ln a^ 

7T 7 



= o 



(42) 



6 



where A = ft log [Hfij / (2Tr)] /(ttMj) measures the 
strength of quantum fluctuations. The master equation 
(131 hence reduces to, 



dt 



+l{p 2 ) 



Mdq 

a 2 



o 2 



M 



W 



(43) 



dp 2 dpdq 

with D(q) = 1 + PXV"(q). As discussed in Refs. [11130], 
the effective diffusion coefficient D(q) should be regarded, 
for thermodynamic consistency, as the first order expan- 
sion of D{q) = (f - pW'tq))- 1 . Equation (43) has been 



derived for arbitrary potentials V(q) in Ref. [3T] (note 
the presence of the incorrect effective potential Veg in the 



latter [H]). The stationary solution of Eq. (43) reads, 

r,2 ' 



W(q,p) = iVexp 



■w* + ^v ( .)'-t 



(44) 



For the special case of a harmonic potential, V(q) = 
Mu! 2 q 2 /2, we can immediately identify the coefficient, 

M 2 lj\ , s 

~kT~ X - ^ 



B = —- 



where we made use of Eq. (24). It is important to real- 



ize that the semiclassical equation ( 43 1 is valid both in 



the high-temperature, kT ^> tvy, and low-temperature, 
kT <C fir/, regimes. I n th e high-temperat ure limit A = 
h 2 /(12MkT) and Eq. (If} reduces to Eq. {27}. 



For general potentials V(q), the deviation from the 
Gibbs state takes the form, 



AH, 



A P - + ^V\q) 2 



(46) 



.4 



(47) 



2 2 

where the coefficient A is given by, 

1 I 

]m~ m 

and is independent of the system. In the semiclassical 
strong-coupling domain, the deviation from the Gibbs 
state can thus be determined for arbitrary potentials and 
not only for the harmonic oscillator. We mention that 
in the high-temperature limit, the semiclassical Hamilto- 
nian of mean force has been examined for a free particle 
and a harmonic dumbbell up to order h 2 in Ref. [35] . We 
note, moreover, that the deviation in position can be di- 
rectly obtained from the simpler quantum Smoluchowski 
equation which follows from the the semiclassical phase- 
space master equation El [2"T1 125] |3T1 132] . 



d_ 
dt 



W(q,t) 



i a 

jMdq 



1 



V'(q,t) + -D(q) 
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W(q,t) , (48) 



where we have introduced the distribution, 

W(q,t) =ps{q,q,t) = J dpW(q,p,t) , (49) 

The corresponding stationary solution is 

A/3 2 



W(q) = ^ 9 exp 



-pV(q) 



V\qf 



(50) 



VI. PHYSICAL INTERPRETATION OF AH S 

A physical interpretation of the difference AH$ be- 
tween the quantum HMF and the system Hamiltonian 
can be given by considering the initial coupling between 
system and bath. By treating the latter process as a ther- 
modynamic transformation, both the change of energy 
and entropy of the quantum system during the coupling 
can be determined [33] . We assume that system and bath 
are initially decoupled, and each in thermal equilibrium 
at temperature T. The total density operator is thus 
given by the direct product, 



exp(-PH s ) exp(-pHB) 

Qy 

Zq Zr 



(51) 



If the system-bath interaction is switched on quasistati- 
cally, the entropy of the system changes by AS and an 
amount Q of heat is exchanged with the heat bath. Both 
can be evaluated by writing the reduced density operator 
of the system Eq. Q in the form, 



PS = exp [~P(Hs + AH S - F)} 



(52) 



where F = —(1/0) InZ* is the free energy of the system. 
The von Neumann entropy of the coupled system, S — 
-ktr s {p s liips}, is then, 



TS = (U-F+ (AH S )) 



(53) 



where U = (Hs) = t?s{psHs} is the internal energy of 
the coupled system. According to the first law, the heat 
exchanged with the bath during the coupling process is 
Q = AU — W, where the work is equal to the free en- 
ergy difference between uncoupled and uncoupled system 
W = AF [12]. We, therefore, have, 



Q = kTAS - (AH S 



(54) 



where we have made use of the fact that, before the cou- 
pling (7 = 0), Uo = kTS — F . The above result, valid 
for any quantum dissipative system, shows that the aver- 
age difference between the HMF and the bare Hamilto- 
nian of the system is just the difference between kTAS 
and the heat exchange Q during the coupling process. 
For vanishingly small coupling, we find Q = kTAS, in 
agreement with standard thermodynamics. 



VII. CONCLUSION 

For a damped quantum system deviations from stan- 
dard thermodynamics occur in the finite coupling regime. 
With the help of the quantum HMF, we have derived an 
exact expression for the deviation from a thermal Gibbs 
state for a damped quantum harmonic oscillator coupled 
to an Ohmic heat bath. We have obtained useful approx- 
imations in the limit of high and low temperature, and of 
weak and strong coupling. In the semiclassical regime, we 



7 



have, moreover, used the quantum Smoluchowski equa- 
tion to derive an approximate formula for the HMF valid 
for any quantum system. Finally, we have established 
a connection between the deviation from a Gibbs state 
and the thermodynamic change of the system that occur 
during the initial coupling process. Our findings empha- 
size the importance of the HMF in the thermodynamic 
analysis of strong damped quantum systems. 

This work was supported by the Emmy Noether Pro- 
gram of the DFG (contract No LU1382/1-1) and the clus- 
ter of excellence Nanosystems Initiative Munich (NIM) . 



Appendix A 



and 



0(1 + a;) = -C+ — x + 0{x 2 ) \x\<l, (A2) 
6 



where C is the Euler constant. On the other hand, the 
series expansions of the hyperbolic cotangent, needed in 
the computation of the coefficients A and B, Eqs. (25) 
and |26l), are El, 



arcoth(y) = — V O 



\v\ > i 



(A3) 



We collect, for convenience, the series expansions of the 
digamma function used in the evaluation of the position 
and momentum quadratures, Eqs. (15) and (16), [34 , 



ip{x) = ln.T 



1 

2x 



I2x 2 



O 



\x\>l (Al) 



and 



arcoth(l + y) 



\y\ < i 



(A4) 
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